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Abstract 

In this paper we propose a novel optimal allocation problem for perennial 
plants. We do not restrict ourselves to consider only the periods favorable for 
photosynthesis, but analyze the whole life period of a perennial plant. This 
provides more information about the strategies of a plant during transitions 
between the favorable and unfavorable seasons. One of predictions of the 
model is that a plant can begin re-establishment of vegetative tissues some 
time before the beginning of the favorable conditions, so as to come into the 
better environmental conditions with a certain amount of already developed 
vegetative tissues. 

In the second part of the paper we extend the model developed in previous 
sections to the case, when a plant controls not only a distribution of allocated 
carbohydrates, but also a size of a seed. We provide sufficient conditions 
under which the optimal strategy for a plant is to produce as much seeds 
as possible. These results can be applied in particular for colonizing species 
and plants living in open environments. The opposite strategy, namely to 
produce the seeds as large as possible, is also considered. 
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1. Introduction 

The pioneering work [l[ gave rise to a new class of mathematical models of 
plants, based on the methods of optimal control theory in which it is assumed 
that a plant can control itself to maximize its fitness to environment which 
is often identified with the mass of seeds produced by a plant during its life. 

The first models have been devoted to the development of annual plants. 
In them it is assumed that a plant consists of a number of compartments, 
at least of vegetative compartment (leaves, roots, stems) and reproductive 
compartment (seeds and auxiliary tissues), but also can include storage and 
defensive tissues. 

The basic model of this type ji[ results in a bang-bang transition from 
the allocation to vegetative tissues to the allocation to seeds. 

This model of an annual plant has been extended in many directions: in 
particular in (if a model with multiple vegetative compartments has been 
analyzed, in works 0|, (EJ a model with additional physiological constraints 
has been considered that results in periods of mixed growth (where both 
vegetative and reproductive parts of a plant grow simultaneously). The op- 
timal allocation strategies in stochastic environments have been investigated 
in particular in Allocation to defensive tissues was encountered, to cite 
a few, in Q, i- 

The overview of the early works in this field is provided in a paper 



For general overview of the resource allocation in plants see books [10| and 
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In contrast to annual plants, to the modeling of optimal phenology of 
perennials less attention was devoted. Usually behavior of perennials is mod- 



eled in the following way [12j, [13j, [14|: the life-time of a perennial plant is 



divided into discrete seasons, which environmental conditions are favorable 
to the photosynthesis. The model of a plant in every season is continuous 
and is treated with the methods used in annual plant models. To model the 
behavior of a plant between seasons (when the weather is unfavorable) some 
simple transition rules are used that show which parts of compartments are 
saved during the season and which are not. 

The solution of such problems is divided into two parts: firstly the model 
on one season is solved using Pontryagin's Maximum Principle (see, e.g. (l5| ) 
and then one seeks a solution of the whole model using dynamic programming 
method. 

Although these models provide quite interesting qualitative results of be- 
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havior of perennial plants, they have an important disadvantage: a subtle 
qualitative behavior of a plant within the season contrasts to the simple jump 
from the end of one season to the beginning of the next one. 

In this paper we are going to propose a continuous-time model of a peren- 
nial plant. This allows us to describe more precisely the dynamics of a plant 
during seasons with unfavorable for photosynthesis environmental conditions 
and to avoid the introduction of additional parameters for description of 
jumps between seasons. 

In particular one of predictions of the model is that the plant begins 
to generate the vegetative tissues not at the time, when the environment 
conditions are favorable for photosynthesis, but some time before, so as to 
come into the suitable period with developed vegetative tissues. 

In the second part of the paper we investigate a trade-off between size 
and number of seeds. A lot of attention is devoted to this topic in the 
scientific literature. The basic model has been proposed in a seminal work 



16| , where it was assumed that the fitness of the plant is equal to the sum of 
the fitnesses of the descendants. Since that this model has been generalized 
in many directions (for a review see 17] ) • In this framework the optimal size 
is sought depending on the properties of the fitness function. This makes 
possible quite general treatment of size-number trade-offs, but the question 
remains, how to formalize the dependency of the fitness on size and number 
of seeds and how to find the properties of the function, that characterizes 
this dependency. 

Our aim is to investigate the trade-offs between number and size of a seed 
in the context of the optimal allocation models. Within this framework the 
fitness is properly formalized, and we can investigate the the optimal size of a 
seed depending on the properties of photosynthetic rate function, and other 
physiological parameters of a plant that are more distinct criteria than the 
abstract fitness. We provide the analysis for the model developed in Section 
|2] of this paper, but the results are valid also for a number of other optimal 
allocation models. 

We prove that, according to our model of a plant, if the photosynthetic 
rate function is concave (that is, if rate of photosynthesis per unit mass 
decays with increase of a size of a plant), then the seeds have to be as small 
as possible. This behavior is usual in particular for the colonizing species 
(see Section [372]) . 

Our model includes also the possibility of choosing of a sprouting time 
by a seed. As a consequence we obtain the results, concerning qualitative 
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behavior of a plant from the dormancy up to senile stage. 

The outline of the article is as follows: in Section 12.11 we introduce the 
model of a perennial plant. In Section 12.21 we provide its analysis, using 
Pontryagin's Maximum Principle. Then, in Section 12.31 we summarize the 
results of the model, provide a general scheme of plant development and 
consider some special cases (annual and monocarpic plants). 

In Section [3] we consider a trade-off between size and number of seeds. 
The results of the paper are discussed in the Section l3~2l Section [4] concludes 
the paper. 

2. Model, in which a plant controls the free carbonates 

2.1. Model description 

Usually it is assumed that all the allocated photosynthate is immediately 
used for construction of the tissues. In the models, taking into account a 
presence of a storage compartment, a plant can also allocate the resources 
from the storage with some maximal rate (depending on the mass of the 
storage) . 

Such a method ignores that a photosynthate is not immediately allocated 
to certain structures, but exists for some time in the free state. We are going 
to take this effect into account and assume that there exists an intermediate 
stage, when the carbohydrates have already been photosynthesized, but have 
not already been permanently allocated to a given structure. 

Let a plant consist of three parts: vegetative compartment, reproductive 
compartment and nonstructural carbohydrates (free glucose, starch etc.). Let 
X\ (t) be the mass of vegetative compartment at time t, 
X2(t) be the mass of reproductive compartment at time t, 
xs(t) be the mass of nonstructural carbohydrates at time t. 

We model dynamics of a plant by the following equations: 



Here shows the rate of photosynthesis of the plant with the veg- 

etative mass x\ in the optimal environmental conditions, and g{x%) - the 
maximal rate of allocation of nonstructural carbohydrates. It is natural to 
assume that / and g are monotonically increasing and /(0) = g(0) = 0. 




(1) 
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The influence of a climate is modeled by three functions: v : [0, T] — > [0, 1] 
and /i, uj : [0, T] — )■ [0, oo). 

z/(t) shows the dependence of the rate of photosynthesis on the climate 
(z/(t) = if at time t no photosynthesis is possible). 

fi(t) is the deconstruction rate of vegetative tissues per unit mass at time t 
u)(t) is the deconstruction rate due to external factors (spoiling, grazing by 
animals etc.) of the storage parts per unit mass at the time t. 

Note that photosynthesized carbohydrates firstly enlarge the mass of non- 
structural carbohydrates. 

We assume that a plant can control the total allocation rate with the 
control v(t) G [0,1], and allocation rate to the vegetative tissues vi(t) G 
[0,f(t)], consequently the allocation rate to reproductive tissues at time t is 
v 2 (t) = v(t)-v 1 (t). 

At the moment i = 0a mass of the seed and all its compartments is given 
a priori (the question of optimal mass of a seed will be considered in the next 
section): 

x l (0)=x° l , z = 1,2,3. (2) 

The seed dormancy is modeled as an ability of a plant to choose the time 
of sprouting t G [0, T]. For simplicity we assume that a seed cannot spoil, 
and it does not use any resources for life-sustaining activities, thus: 

x i (t o ) = x i (0) = x° i , i = 1,2,3. (3) 

As the aim of a plant we choose maximization of the total yield of seeds 
over the period: 

X2(T) — > max . (4) 

A plant can achieve this aim choosing the appropriate sprouting time to 
and controls v and i>i defined on [t , T] . 

We assume that the functions in the right hand side of equations ([1]) are 
smooth enough to guarantee existence and uniqueness of solutions of ([T|). 
Also we assume that the system ([T]) is forward complete, that is, for all 
initial conditions and all admissible controls the solutions of ([T]) exists for all 
time. From the biological viewpoint it means that it is impossible to achieve 
endless yield in finite time. This assures that the solution of the problem ([T|), 
(J2D, flU) exists. 
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2.2. Model analysis 

For analysis we exploit Pontryagin's Maximum Principle (see, e.g. 
The Hamiltonian of the system (CQ) is defined by: 

H = pi (vi(t)g(x 3 ) - fi(t)xi) + p 2 (v(t) - v ± (t)) g{x 3 )+ , , 

p 3 {v{t)f(x x ) - v(t)g{x 3 ) - u(t)x 3 ) . {0) 

Functions Pi,P2,P3 are so-called adjoint functions. The equations determin- 
ing their dynamics will be given later. We rewrite expression (J5J) in a more 
suitable form 

H = p 3 u(t)f(xx) - Pifi(t)xx + g(x 3 )(v 1 (t)(p 1 - p 2 )+ 

V(t)(p 2 -Ps)) ~p 3 U){t)x 3 . 

For computing of the boundary values of functions pi, i = 1, ... ,3 the 
following function is needed 

3 

/(to, x(t ),x(T)) = -X x 2 {T) + W*i(to) ~ 4) + ^o- (7) 

i=i 

Here A > 0, A 4 < 0, \ G R, i = 1, 2, 3. 
Equations for the adjoint function p are as follows 

px = pxfi(t) - p 3 v(t)§L(xx), 

P2 = 0, (8) 

^3 = -^(x 3 )(Vx(px ~P2) + V(p 2 -p 3 )) +PsU(t). 

The corresponding boundary conditions are 

Px(T) = 0, 

P2(T) = A , (9) 
Ps(T) = 0. 

Initial conditions are 

dl 

= TwTT = A - * = 1.2,3. (10) 

OXi{t ) 

Possibility to choose optimal to provides us with one more equation 

H(to) = ~ = -A 4 . (11) 
6 



Also the condition of complementary slackness must hold 

A 4 • t = 0. (12) 
To obtain the values of v,v±, we solve the problem 
H — > max, < v < 1, < V\ < v. 

It is not hard to check that its solution is given by 

1. If P2 — p-3 > 0, then v — 1, and 

{v if p 1 - p 2 > 0, 

ifp!-p 2 <0, (13) 

e [0, v] if pi-p 2 = 0. 

2. If p 2 — p-3 = 0, then 

Pi — P2 > =^ W = 1 , V\ — V 

Pi — P2 = =^ v, Vi - every possible controls (14) 
Pi — P2 < =>• f - every possible control, v x = 

3. If P2 — Pz < 0, then 

• if £>i — P2 < then ?; = V\ — 0. 

• if pi — p 2 > then 

Pi — p 3 < =^ u = Vi = 

Pi ~ P3 — v - every possible control, V\ = v (15) 

Pl — P3 > =^ V — V\ — \ 

We introduce three main periods characterized by different values of con- 
trols: 

1. Vegetative period: pi > max{p 2 ,P3}- In this case v — Vi = 1, that is 
the vegetative parts are being constructed with maximal rate. 

2. Reproductive period: p 2 > ma,x{pi,p 3 }. In this case v = 1, v i = and 
reproductive tissues are being constructed with maximal rate. 

3. Storage period: ps > max{pi,p 2 }- In this case v — v± — and new 
tissues are not being constructed. 
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We are going to analyze these periods more deeply and find out in what 
order these periods can arise in a life of a plant. To this end we investigate 
equations (jHJ) from the end of the life of a plant. 

If A = 0, then from and flS} we obtain that pi = on [to,^], from 
which it follows that all the controls are possible. 

Let A > 0. We can take in this case A = 1. 

Controls v and v 1 maximize the value of (vi(pi—p 2 )+v(p 2 —p 3 )), therefore 
for optimal v,v± it holds that 

(vi(p 1 -p 2 )+v(p2-p3))>0. (16) 

Note that in case, when oo(t) = (that is, if storage parts cannot be 
destructed due to external factors) this inequality and monotonicity of g 
imply that p 3 is an non-increasing function on [to ; ^]- 

Let us analyze the behavior of Lagrange multipliers pi and values of con- 
trols at the neighborhood of the time T. 

From conditions we see that at the moment t — T it holds P2{T) > 
max{pi(T),p 3 (T)} and from continuity of pi, i — 1,2,3 it follows that the 
last period of the plant life is always a reproductive period. 

From equations (fTBl) and dHJ using monotonicity of g and inequality uj > 
we have that if for some r G [t , T] p 3 {j) < 0, then p 3 (t) < for all t G [r, T], 
which contradicts to Thus, p 3 > on [t , T]. Analogously one can prove 
that pi > on [t ,T]. 

During a reproduction period equations (jSJ) after substitution of corre- 
sponding controls take a form 

h=p 1 (t)»(t)- P3 (t)v(t)£ i {x 1 (t)), 

h = 0, (17) 

P3 = -^- 3 (x 3 (t))(p 2 (t)- P3 (t))+Uj(t)x 3 (t). 

From these equations it follows that 

sup{r G [0, T] : p x {r) = p 2 {r)} > sup{r G [0, T] : p 3 {r) = p 2 (r)}. 

This tells us that before reproduction period the storage period is impossible. 
Let ^ be the last (nearest to T) moment, such that Pi{t\) = p 2 (ti). 

If the climate conditions (functions fi and v) are such that Pi(t) = p 2 {t) 
for all t G [t s , ti] for some t s < t\, then according to ( [TBI plant can have the 
period with mixed control v\ G [0,v] for t G [t s ,ii]. Although this possibility 
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cannot be excluded in general, but such mixed controls can arise only due to 
very specific climate conditions and we do not separate it as a special period 
of plant life. 

If pi is increasing at t — tx, then one can distinguish one more repro- 
ductive period [ti — s,ti) for some s > 0. Throughout this paper we follow 
the agreement to combine all such periods together with stages with mixed 
controls between these periods into one reproductive period. 

Let pi be decreasing. Then for some time interval preceding to the re- 
productive period we have p\ > P2 > p 3 and therefore on this time interval 
a plant has a vegetative period. We will call it period 1.1, in contrast to 
1.2 characterized by relation pi > p 3 > p 2 (this distinction is useful for 
monocarpic plants). 

Equations ()8]) in the vegetative period take form: 



There are 2 possibilities for the plant behavior before period 1: either it 
will have one more 0-th period (if pi decreases lower than p 2 while it remains 
true that p 3 < p 2 ), or it will exist t 2 < t\\ p 3 {t 2 ) = P2(t 2 ). As mentioned 
before, we neglect the possibility of mixed controls and consider the case 



In this case period 1.2 characterized by p\ > p 3 > p 2 precedes the period 



To understand the difference between periods 1.2 and 1.1, let us consider 
the case, when u = (which implies, as was mentioned earlier, that p 3 is 
non- increasing) . This implies that before period 1.2 the reproduction periods 
are not possible (p 3 > p 2 ) and consequently the plant exploits monocarpic 
strategy. 

In more general case, when u ^ both periods 2 and 3 can precede the 
first period, or all the previous life of a plant can consist of one vegetative 
period. In the first plant possesses one more reproduction period, 

which has been already analyzed. If before vegetative period there is no 
other period, then the plant is annual. 

Let now the third period precedes to the first one. Then there exist t 3 : 
ti < t 3 < t 2 , such that pi increases on [t^,t 3 ] (due to the unfavorable climate 
conditions) and ^1(^4) = ^3(^4). We separate period between £4 and £3 in the 



p 1 = p 1 (t)fi(t)-p 3 (t)v(t)£ I (x 1 (t)), 

P2 = 0, 

^3 = -^-(X 3 {t)){p 1 (t)- P3 {t))+U(t)x 3 . 



(18) 



1.1. 
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season 1.2.1 (pi > p 3 > p 2 , but pi decreasing), which distinctive feature is 
that although the climate conditions are not comfortable for photosynthesis 
a plant anyway allocates some part of stored resources to the construction of 
the vegetative tissues, so as to come into the better conditions with certain 
amount of already developed vegetative mass. 

Now let there exist some r: pi(t) < p 3 {t) for all t G [r, t 4 ). Then a plant 
enters a storage period. The corresponding equations ([H]) take a form 

Pi = Pi(t)^(t) - p 3 (t)v(t)£L( Xl (t)), 
P2 = 0, (19) 
p 3 = u(t)x 3 . 

If the climate conditions are unfavorable for all t < £ 4 , that is, Pi(t) < 
p 3 {t) for all t G [0,£ 4 ), then the first period of time is only the storage of 
allocated photosynthate (this is hardly possible because a seed has a possi- 
bility to stay this period in dormancy). If it is not the case, then there exist 
some moments t e , t 5 , t 6 < £ 5 < t 4 , such that p\ is decreasing on [t e ,t 5 ] and 

Pi(k) = Pz(h)- 

We separate the period (t 5 , i 4 ), which we call period 3.1 (when the climate 
conditions are disadvantageous and all the allocated material is stored), and 
time-span (^^5) called period 3.2 (when the climate conditions are kindly, 
but all the allocated material is anyway stored "preparation to the unfavor- 
able climate conditions"). 

Both reproductive and vegetative periods can precede to the storage pe- 
riod. It depends on the climate conditions and values of x°. More precisely, 
from (HI]), (ED and (HO]) we have: 

H(t ) = g(x 3 (to)) M£ )(Ai - A 2 ) + v(t )(X 2 - A 3 )) + (2Q) 
A 3 ^(*o)/(^i(*o)) - Ai/x(t )zi(£o) + \ 3 u(t )x 3 (t ). 

If to > 0, then according to the condition of complementary slackness ( Tl"2"D 
H(t ) = 0. 

It seems that in general we cannot say more about the time of propagation 
and type of the first period. The reason is that one can choose the values of 
the initial parameters that are biologically inadequate and consequently ob- 
tain unrealistic predictions. For example, if the climate conditions are chosen 
to be unfavorable for photosynthesis throughout all the time-interval [0,T], 
then the model is inapplicable, because the strategy to stay in dormancy all 
the period is not allowed in the model. 
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Figure 1: Stages of perennial plant development 



To exclude such biologically irrelevant behavior, we consider in the next 
subsection the case, when the first period after sprouting is vegetative period. 

2. 3. Predictions of the model 

In this section we conclude our investigations. The development of the 
plant according to the model ([I]) consists of 3 main periods that can be further 
subdivided into subperiods and that can follow each other as depicted in the 
Figure [TJ 

These periods are: 

• - Dormancy. 

• 1 - Vegetative period. 

• 3.2 - Preparing for the unfavorable climate conditions. 

• 3.1 - Life in unfavorable climate conditions. 

• 1.2.1 - Allocation to vegetative tissues as a preparation for the climate 
conditions favorable for the photosynthesis. 

• 1.2.2 - A vegetative period that is important for the characterization 
of monocarpic plants (see below). 

• 1.1 Allocation to vegetative tissues before reproduction. 

• 2 - Reproduction. 

Important special cases of this scheme are: 

1. Annual plant with multiple reproduction periods: see Figure [3] (period 
1 we identify with 1.1 for simplicity). Multiple reproduction periods 
appear because of losses of vegetative mass due to external factors that 
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Figure 2: Life-stages of annual plants 
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Figure 3: Life-stages of monocaprics 



are modeled by the function fj,. This particular case has been analyzed 



in the early work [181 ] . If fi = 0, then the multiple reproductive periods 
for annual plants are not possible. 
2. Monocarpic plants. Transition 1.1 — > 1.2.1 is not possible, see Figure 
I3J Sufficient (but not necessary) condition for a plant to be monocarpic 
is negligibility of u (in particular, if uj = 0), in other words, the mass 
of storage cannot decrease due to the external factors. 



3. Optimization of the seed mass 

In the previous sections we have defined fitness of a plant as a mass of 
reproductive tissues produced by the plant during its life. To maximize the 
fitness, a plant controls the allocation of photosynthate. 

However, it is well-known, that the fitness depends (for plants that prop- 
agate themselves with the help of the seeds) crucially on the quantity (and 
size) of seeds, that a plant produces. Current models of optimal allocation do 
not provide this information, and a mass of a seed is treated as an external 
parameter. 

3.1. Modeling of seed mass optimization and main result 

In this section we consider the model of optimal allocation developed in 
previous sections and consider the optimal seed size problem for it. 

Let yo = (ui, y®, y®) be the total mass of the seeds (the vector consisting 
of the masses of three components of a plant) that has to be divided between 
a seeds, a £ [1, oo) and a can be either natural or real number. 
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We assume that mass of each seed is s = — . 

a 

The equations determining dynamics of a plant are as follows 

±1 = vx(t)g(x 3 ) - n(t)xi, 

X2 = (v(t)-v 1 (t))g(x 3 ), / 21 x 
^3 = v(t)f{xi) - v(t)g(x 3 ) - u(t)x 3 , 
x{0) = >■ 

Here x(0) = (xi(0), x 2 (0), x 3 (0)). 

The aim of a plant we define as maximization of the total mass of repro- 
ductive tissues produced by all direct descendants: 

max Q a = ax n (T). (22) 

0<v(t)<l, 0<vi(t)<v(t), o6[l,oo) 

In contrast to the problem (JTJ), (J4]) with the fixed mass of the seed, the 
problem (T2~Tj) . ( )22|) may have no solution (optimality might be achieved for 
a — > 00), that is, no admissible controls v, v\ and parameter a generate the 
optimal value of ax n {T). 

Recall that a function / is called concave on the set M, if Vx, y G M, 
V« G [0, 1] it holds inequality 

f(ax + (1 - a)y) > af(x) + (1 - a)f(y). (23) 

If the inequality ( 1231) holds with < instead of >, then the function / is called 
convex on the set M. 

Important is the case, when / and g are concave functions, that is, the 
rate of photosynthesis and maximal speed of chemical reactions in a plant 
are saturated with the growth of the mass of a plant (due to self-shading of 
leaves, nutrient depletion in the soil etc.). 

We have the following result: 

Proposition 1. Let f,g : R + — > R + be concave, /(0) = g(0) = 0, a G [l,oo). 
Then Q a = ax2(T) increases, when a increases (x2(T) is a solution of a 
problem 

Proof. The problem ( 12T|) .( 122|) can be written in equivalent form, using new 
variables yi{t) := axi(t), i = 1,2,3. Then we have: 

2/1 = Vi{t)ag{^) - fi(t)yi, 
y 2 = (v(t)-v 1 (t))ag(yf), 

h = v{t)af(K) - v(t)ag(f) - u(t)y 3 , ^ 
2/(0) = y - 
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The corresponding maximum problem is: 

max Q = y 2 (T). (25) 

0<v(t)<l, 0<vi(t)<v(t), aG[l,oo) 

Now the problem is similar to (CD), (jlj), but with a/(— ) and ag(—) instead 
of f(x\) and g(x 3 ). 

Using concavity we have: f{f) = f{\y x + 2=1 . 0) > ±/( yi ) + 2=1/(0) = 

Thus, for every yi(t) > 0, a > 1 it holds afip^-) > f(yi(t)) and therefore 

a/(^^) and agi^p-) are nondecreasing in a and sup ae [ 1)00 ) af(^jp-) and 

sup ae[ i j0o) ag(^) yields, when a ->■ oo. 

Define the value of 2/2 (r) for a given a as y2(T,a). Now take arbitrary 
n > a and consider a system 

& = ui(t)a#(^) - //(t)2/i, 

k = (v(f)-«i(*)M?), ,otf 
2/3 = v{t)nf{%) - v(t)ag(f) - u(t)y 3 , W 

2/(0) = 2/0- 

The solution of this system subject to optimality condition (124]) we denote 
y(t). If i/(0) > 0, then from n/(^) > of(^) we have that £ 3 (0) > y 3 (0, a) 
and therefore there exists t* > : ^(t) > y 3 (t,a) V£ G [0,f). Hence 
&(*) > V3(t,a) and a#(^) > ac/(^^) for t G (0,**). Let i> and ui 
are optimal controls for a system ( |24"1) . There exist controls < v < v, 
< Vi < Vi for a system ( 126]) . such that v(t)ag(^^-) = v(t)ag( y3 ^' a ' ) and 

ik{t)ag{&&) = vy^agi^). Consequently, Vi (t) = yi (t,a), t G [0,t*), 
z = l,2. Constructing v, v\ for all t G [0,T], we obtain that 2/(T) = y(T,a) 
and thus the optimal trajectory of a system ( 126!) produces no less seeds than 
the best trajectory of ( 1241) . 

Analogously, the output of the following system is not less than that of 
the system (126]) : 

jfi = v^ngi^) - (i{t)yi, 

V2 = (v(t)-v 1 (t))ng(%), , - 

m = v{t)nf{%) - v{t)ng{^) - u(t)y 3 , [A<) 
2/(0) = 2/0- 

Hence y2(T,a) is nondecreasing in a. □ 
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It follows from the previous proposition that in the case, when / and 
g are concave the best strategy for a plant is to produce as much seeds as 
possible. 

Similar argument shows that for convex functions /, g the optimal mass 
of the seed has to be as large as possible (without additional restrictions on 
quantity of seeds a = 1). 

Remark 1. Assuming that f and g are concave and continuously differ- 
entiable in the neighborhood of one can investigate \im a ^ 00 y 2 (T,a) using 
formulas 

) -> /' (0)|/i (t) and ag(^-) g'(0)y 3 (t), when a -+ oo. 
a a 

The solution of the obtained linearized problem (with corresponding optimality 
condition) provides the '"theoretical" ' upper bound for the fitness of the plant 
under consideration. 

Remark 2. Note that if both f and g are linear, then from (Efy it follows 
that the yield of a plant does not depend on the mass of the seeds. 

The results of this section can be analogously stated also for many other 
optimal allocation problems. In particular they hold for the model 



Xl{t) = 


(1 


x 2 (t) = 


u(i 


xi(0) = 




x 2 (0) = 


o, 



• (28) 

r 2 (0) = 0, 
with optimality condition 

max Q = x 2 (T). (29) 

0<w(*)<l 

Here Xi(t) and x 2 (t) are the masses of vegetative and reproductive tissues 
at time t respectively and / is a monotonically increasing function. 

3.2. Discussion of results 

We see that the optimal size of a seed depends crucially on the form of 
the functions / and g. For concave functions that are often used to take 
into account the self-shading, boundedness of resources etc. (see e.g. [I3|) 
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we have proved that according to our model seeds have to be as small as 
possible. 

For the plants living in open environments and for species occupying early 
phases in succession (colonizing species) the assumption of concavity is not 
an oversimplification. The behavior that our model predicts, namely that 
the optimal strategy is to produce a vast amount of small seeds is usual for 



these species [19 



However, in the closed and shady environments, under mineral shortage, 
or if there is a strong competition with the established vegetation, the rate 
of photosynthesis per unit mass can increase with the increasing of a mass of 
a plant, that is, the function / is convex on some [0,p], p > and the seeds 
cannot be too small. These predictions are, in general, in accordance with 



experiments [20fl, [2l|, but see [22 



Note that similar optimal allocation models can be used also to study 



life-strategies of the animals, see 23|. In this case it is worth mentioning 



that there exist species, for which the used in Section [3] definition of fitness 
(the fitness of a parent is equal to the sum of the fitnesses of all descendants) 
is not realistic. For example, in recent paper [24| the experimental studies on 
keelback snakes (Tropidonophis mairii, Colubridae) were provided according 
to which an increase in clutch size increases the offspring size due to the 
fact that the uptake of water by the snake egg depends upon the number 
of adjacent eggs. For this case another models have to be constructed and 
analyzed. 



4. Conclusion 

In the present paper we have developed the model that describes the 
optimal allocation strategies of a perennial (as well as annual) plant during 
all the stages of its life. The model was analyzed with the help of Pontryagin's 
Maximum Principle, and as a consequence the division of a life of perennial 
plant (Figure [1]) was obtained. The models of monocarpic perennial plant as 
well as the classic models of annual plant with one or several reproduction 
periods have been obtained as a special case of the general model. 

In the second part of the paper we analyzed the trade-off between size 
and number of seeds for the optimal allocation problems. We have provided 
the sufficient conditions that ensure that an optimal strategy is to produce 
as much (or as less) seeds as possible. The applicability of results has been 
discussed in the Section [3T2l 
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